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Abstract 

A very famous result of gauge/gravity duality is the universality of the ratio of 
shear viscosity to entropy density in every field theory holographically dual to 
classical, two-derivative (Einstein) gravity. We present a way to obtain deviation 
from this universality by breaking the rotational symmetry spontaneously. In 
anisotropic fluids additional shear modes exist and their corresponding shear 
viscosities may be non-universal. We confirm this by explicitly calculating the 
shear viscosities in a transversely isotropic background, a p-wave superfluid, 
and study its critical behavior. This is a first decisive step towards further 
applications of gauge/gravity duality to physical systems. 
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Gauge/gravity duality [1] provides a novel method for studying strongly- 
coupled systems at finite temperature and densities. As such it is expected 
to have useful applications to the quark-gluon plasma as well as to condensed 
matter physics. The duality also allows to calculate physical observables in the 
real-time formalism, e. g. retarded correlation functions [2, 3] and hydrodynamic 
transport coefficients [4]. The most important result of these studies is the 
universality of the ratio of the shear viscosity -q to the entropy density s in every 
field theory holographically dual to classical, two-derivative (Einstein) gravity 
theory [5, 6, 7, 8]. The universal number rj/s = 1/A-k (in natural units) fits 
surprisingly well the measured quantities in strongly-coupled real- world systems 
such as the quark-gluon plasma and cold fermions at unitarity, see e. g. [9]. The 
classical, two-derivative gravity theory is mapped to a gauge theory with large 
rank of the gauge group (large N c ) and to a large 't Hooft coupling A. 

Since the ratio of shear viscosity to entropy density usually depends on the 
temperature and is non-universal in real-world systems, it is important to find 
deviations from the universal behavior in order to make further contact with 
real-world systems. In this letter, we present a deviation from the universal 
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ratio at leading order in N c and A. So far a deviation from the universality only 
occurs in subleading corrections of N c and/or A, see e.g. [10, 11]. 

The universality of 77/s strongly depends on the shear mode transforming as 
a helicity two state under the rotational symmetry. Since it is the only helicity 
two state in the system, the shear mode decouples from the other modes and 
it can effectively be written as a minimally coupled scalar. The universality of 
77/s is then related to the universality of the coupling of this minimally cou- 
pled scalar. The coupling is always given by the gravitational coupling constant 
since the shear mode is a graviton [8] . We circumvent this important assumption 
that the shear mode is a helicity two state by breaking the rotational symme- 
try spontaneously. A fluid with spontaneously broken rotational symmetry is 
anisotropic and more than one shear mode exists. Since some of these modes 
do not transform as helicity two states, the corresponding shear viscosities can 
be non-universal [12, 13]. 

The viscosity, including shear and bulk viscosity, in anisotropic fluids is de- 
scribed by a rank four tensor with in general 21 independent coefficients [14]. In 
the following we study systems with p-wave symmetry (transversely isotropic). 
This reduces the independent coefficients in the viscosity tensor to five. Two of 
these independent coefficients are shear viscosities [15]. We explicitly calculate 
these two shear viscosities by using the recipe for the holographic calculation of 
the retarded correlators [2, 3] and Kubo formulae [16] which are given by 

r/i = — lim -ImGfj(w,0) with i e {xy,xz,yz} , (1) 

where the x-axis is the preferred direction and rj xy (= T] xz ), rj yz are the two 
independent shear viscosities. The retarded correlation function of the energy- 
momentum tensor T^ v is defined by 

G& w (w,0) = -iJ dt&xe wt 6(t){[T l0 {t,x),T kl (0M) ■ (2) 

We present more details on anisotropic shear viscosities in the appendix below. 

In the holographic context, the spontaneous breaking of symmetries by black 
holes developing hair was first achieved in [17] and later used to construct 
holographic superconductors/supcrfluids by breaking an abelian symmetry [18]. 
Along this line also p-wave superconductors/superfluids have been constructed 
[19] and gave rise to the first string theory embeddings of holographic super- 
conductors/superfluids [20, 21, 22]. In p-wave superconductors/superfluids, also 
the spatial rotational symmetry is spontaneously broken in addition to the inter- 
nal abelian symmetry. In order to obtain the effects of spontaneous rotational 
symmetry breaking for the energy-momentum tensor, which determines the hy- 
drodynamics of the systems, we have to take the back-reaction of the superfluid 
into account. This was obtained e.g. in [23]. On the gravity side, the p-wave 
superfluid state corresponds to an asymptotically AdS black hole which carries 
vector hair. 

In order to construct p-wave superfluid states, we consider SU (2) Einstein- 
Yang-Mills theory in (4 + 1) -dimensional asymptotically AdS space as in [23]. 
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The action is 



a 2 



{R-K)-^F; V F^ V 



(3) 



where K5 is the five-dimensional gravitational constant, A = — -p- is the cosmo- 
logical constant, with L being the AdS radius, and a = K§/g the ratio of the 
gravitational constant k 5 to the Yang-Mills coupling constant g. The SU(2) 
field strength F* v is 

F^ = d^A a u -d v A; + e abc AlAl, (4) 

where fi, v = {t,r,x,y, z}, with r being the AdS radial coordinate, and e abc is 
the totally antisymmetric tensor with e 123 = +1. The A® are the components 
of the matrix- valued gauge field, A = A^T a dx^, where the r a are the SU{2) 
generators, which are related to the Pauli matrices by r a = a a /2i. 

Following ref. [19], to construct charged black hole solutions with vector hair 
we choose a gauge field ansatz 

A = 4>(r)T 3 dt + w(r)r 1 dx . (5) 

The motivation for this ansatz is as follows. In the field theory we introduce 
a chemical potential for the U(l) symmetry generated by r 3 . We denote this 
U(l) as (7(1)3- The gravity field dual to the U(l) 3 density is Af, hence we 
include A^(r) = tfi(r) in our ansatz. We allow for states with a nonzero (Jf), 
so in addition we introduce A].{r) = w(r). With this ansatz for the gauge 
field, the Yang-Mills stress-energy tensor is diagonal. Solutions with nonzero 
w(r) preserve only an SO(2) subgroup of the 5*0(3) rotational symmetry, so 
our metric ansatz respects only 50(2). Furthermore, given that the Yang-Mills 
stress-energy tensor is diagonal, a diagonal metric is consistent. Our metric 
ansatz is 

ds 2 = - N{r)a(r) 2 dt 2 + ^ydr 2 + r 2 /(r)- 4 d.x 2 
+ r 2 /(r) 2 (dy 2 +dz 2 ) , 

with N(r) = -^P + £. For our black hole solutions we denote the position 
of the horizon as tv The AdS boundary is at r — > 00. 

Inserting our ansatz into the Einstein and Yang-Mills equations yields five 
equations of motion for m(r), cr(r), f(r), 4>{r), w{r) and one constraint equation 
from the rr component of the Einstein equations. The dynamical equations can 
be found in [23]. Using scale transformations, we can set the boundary values 
of both cr(r) and f(r) to one, so that the metric will be asymptotically AdS. 

A known analytic solution of the equations of motion is an asymptotically 
AdS Reissner-Nordstrom black hole, which has <p{r) = /j, — q/r 2 , w(r) = 0, 

°(r) = f(r) - 1, and N(r) = (r 2 - ^ + ^#), where TOo = ^ + ^# and 

q = fir 2 ,. Here \i is the value of <j>{r) at the boundary, which in CFT terms is the 
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U(l)3 chemical potential. Since w = 0, i.e. (Jf) = 0, this solution corresponds 
to the normal phase of the system. To find solutions with nonzero w(r) we 
resort to numerics. We solve the equations of motion using a shooting method 
(see [23] for details). 

We calculate the shear viscosity in this background via the Kubo formulae 
by considering two point functions of the energy-momentum tensor and the 
currents 1 . As described in [2, 3, 24], these two point functions are determined 
by considering fluctuations of the metric h^ u and the gauge fields about this 
background. For these fluctuations we choose the ansatz 

Taking the gauge freedom into account, we may set the modes a% and h^ r for 
jU € {t, x, r} to zero. This leads to eight constraints in addition to the equations 
of motion for the 22 dynamical fields. This reduces the number of the physical 
modes of the system to 14. 

The shear viscosities are determined by the zero momentum correlation func- 
tions. Thus we only consider time dependent fluctuations, i.e. k — 0. These 
fluctuations can be characterized by their transformation properties under the 
unbroken SO (2) rotational symmetry: 

hclicity 2 : h yz , h yy — h zz 
helicity 1 : h yt , h xy , h yr ; a a y 

hzy> hxzi h zri d z (8) 

helicity : htt, h yy + h zz , h xx , h xt , h xr , h tr , h rr ; 
a t ,a x ,a r . 

The 14 physical modes decouple in several blocks. The first block contains 
the usual two physical helicity two modes including the shear mode h yz . In 
the helicity one block there are eight physical modes which can be split in two 
blocks by the residual 5*0(2) rotational symmetry. In this block the additional 
shear modes h xy and h xz appear. The helicity zero block contains additional 
four physical modes. 

Since the first shear mode h yz is the usual helicity two mode, the general 
proofs of universality [5, 6, 7, 8] apply and we obtain r\ yz js = 1/Air. The second 
shear mode h xy is a helicity one mode and couples to other physical helicity 
one modes. Thus the assumptions of the proofs of universality [5, 6, 7, 8] are 
not satisfied and the shear viscosity can be non-universal as we will see in the 
following. 



lr The two point functions of the currents have to be considered in addition to the two point 
functions of the energy-momentum tensor since there are interactions between the energy- 
momentum tensor and the currents. 
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Due to the residual Z2 symmetries, the four physical helicity one modes de- 
couple into two blocks of coupled differential equations in the case of k = 0. 
The first block contains one physical mode which is determined by the dynam- 
ical fields h yt and a y and the constraint h yr — 0. This mode determines the 
electrical [25] and thermal conductivity and the thermoelectric coefficients. The 
three coefficients are related by Ward identities (see e. g. [26]). The second block 
contains the three physical modes a y , a 2 , and h xy . We are interested in the sec- 
ond block which contains the shear mode ^ — h xy /(r 2 f 2 ) = h x . The linearized 
equations of motion for the modes in the second block are given by 



,,. ,'1 4r 6/' 

°=* +i r + N + T 



ra 2 4>' 2 
3No 2 
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2a 2 

r 2 J2 



ibjwcpa; 
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+ < 



f 4 w 2 uj 2 

" r 2 N + N 2 a 2 



+ 



N 2 o 2 



+ 



iu>4> 
N 2 o 2 



where the prime denotes the derivative with respect to the radial variable 
r. We solve the coupled differential equations numerically and determine the 
retarded correlator G xy xy using the recipe of [2, 3, 24], 



G xy . xy (^,0) - +^2 r ' 



— counter terms . 



(10) 



Using the Kubo formulae (f ), we numerically determine the ratio of shear viscos- 
ity r) xy to entropy density s, which is given by the Bekenstein-Hawking entropy 
of the black hole s = 2irrf l /n 2 .. 

In fig 1 we compare our numerical results for the ratio of the shear viscosity 
rj xy to the entropy density s with the universal behavior of the shear viscosity 
rj yz for different values of the ratio of the gravitational coupling to the Yang- 
Mills coupling, denoted by a. We see that in the normal phase T > T c , the 
two shear viscosities coincide as required in an isotropic fluid. In addition, 
the ratio of shear viscosity to entropy density is universal. In the supcrfluid 
phase T < T c , the two shear viscosities deviate from each other and r\ xy is 
non-universal. However it is exciting that the KSS bound on the ratio of shear 
viscosity to entropy density [5] is still valid. 

The difference between the two viscosities in the superfluid phase is con- 
trolled by a. In the probe limit where a = 0, the shear viscosities also coincide 
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in the supcrfluid phase. By increasing the back-reaction of the gauge fields, 
i. e. rising a, the deviation between the shear viscosities becomes bigger in the 
superfluid phase as shown in fig 1. If a is larger than the critical value found 
in [23] where the phase transition to the superfluid phase becomes first order, 
the shear viscosities are also multivalued close to the phase transition as seen 
in fig. 2. Since there is a maximal a denoted by a max for which the superfluid 
phase exists, we expect that the deviation of the shear viscosity r\ xy from its 
universal value is maximal for this a max . Unfortunately numerical calculations 
for large values of a are very challenging such that we cannot present satisfying 
numerical data for this region. It is interesting that also the deviations due to 
A and N c corrections are bounded. In this case the bound is determined by 
causality [27]. 

For a smaller than the critical value where the phase transition is second 
order, we may study the critical behavior of the ratio of the shear viscosities 
to entropy density close to the phase transition. Due to universality, r\ yz js is 
constant and does not change on both sides of the phase transition, while r) xy /s 
is only constant in the normal phase, but has a different critical behavior in the 
superfluid phase. Let us consider the critical exponent related to rj xy /s and its 
dependence on a. From our numerical data we obtain the critical behavior 



in the superfluid phase close to the phase transition. It is interesting that the 
critical exponent /? does not change with a. A more precise statement can be 
made if it is possible to determine this critical exponent analytically. We suggest 
that this can be achieved by an expansion in the order parameter (Jf ) and in 
the back-reaction controlled by a. So far, the expansion is known in the probe 
limit a = only [28]. However we expect that this expansion may be extended 
to small values of a. We plan to study the critical behavior of this system in 
more detail in the future and to make contact with the theory of dynamical 
critical phenomena [29]. 

To our knowledge this is the first non-universal behavior of the shear viscosity 
to entropy density ratio in a theory which is holographically dual to classical, 
two derivative Einstein gravity theory. This non-universal behavior is due to 
the fact that the corresponding gravity fluctuation h xy transforms as a helicity 
one state and therefore couples to the fluctuations of the gauge field a*' 2 . We 
expect that similar results can be obtained for every gravity background which 
spontaneously breaks the rotational symmetry. 

Appendix A. General Remarks on viscosity in anisotropic fluids 

Here we explain why the Kubo formulae (1) are appropriate for calculating 
the shear viscosities r) xy (= r\ xz ) and rj yz even in a p-wave superfluid. 

In general, viscosity refers to the dissipation of energy due to any internal 
motion [30]. For an internal motion which describes a general translation or a 




with /3 = 1.00±3% 
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Figure 1: Ratio of shear viscosities r\ yz and rj xy to entropy density s over the reduced tem- 
perature T/T c for different values of the ratio of the gravitational coupling constant to the 
Yang-Mills coupling constant a: The color coding is as follows: In yellow, r] yz /s for all values 
of o; while the curve for r) xy / s is plotted in green for a = 0.032, red for a = 0.224 and blue for 
a = 0.316. The shear viscosities coincide and are universal in the normal phase T > T c . How- 
ever in the superfiuid phase T < T c , the shear viscosity r\ yz has the usual universal behavior 
while the shear viscosity rj xy is non-universal. 



general rotation, the dissipation is zero. Thus the dissipation depends on the 
gradient of the velocities u M only in the combination u^ v = | {d^u v + d^u^), 
and we may define a general dissipation function S = \rf ,v p u^, v Ux p , where 
rf vXp defines the viscosity tensor [14]. Its symmetries are given by 

The part of the stress tensor which is dissipative due to viscosity is defined by 

TZ s = -^ = -ir Xp ux P . (A.2) 

We consider a fluid in the rest frame of the normal fluid u* = 1 . To satisfy the 
condition of the Landau frame WmTj*. = 0, the stress-energy tensor and thus the 
viscosity has non-zero components only in the spatial directions i,j = {x,y,z}. 
In general only 21 independent components of rjijki appear in the expressions 
above. 

For an isotropic fluid, there are only two independent components which are 
usually parameterized by the shear viscosity rj and the bulk viscosity £. The 
dissipative part of the stress tensor becomes T^ lss = —2r/(uij — ^8iju\) — C,u\8ij 
which is the well-known result. 

In a transversely isotropic fluid, there are five independent components of 
the tensor rjiju- For concreteness we choose the symmetry axis to be along the 
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Figure 2: Ratio of shear viscosities r\ yz (blue) and r] xy (red) to entropy density s over the 
reduced temperature T/T c for a = 0.447, which is bigger than the critical value where the 
phase transition becomes first order: The shear viscosities coincide in the normal phase T > T c 
and are universal. In the superfluid phase r\ xy is non-universal. Close to the phase transition, 
it is multivalued as expected for a first order phase transition. 



x-axis. The non-zero components arc given by 



Vyyy 



Cy 



7]xxyy — TJxxzz 
Vyzyz = Vyz j 



A , Vyyzz — Cy 



Vyz 



(A.3) 



Tjxyxy — Vxzxz — V%y 



The non-zero off-diagonal components of the stress tensor are given by 



iTjxyUxy i 
'^^lyz^yz • 



T 



diss 



(A.4) 



So far, we only considered the contribution to the stress tensor due to the 
dissipation via viscosity and found the terms in the constitutive equation which 
contain the velocity of the normal fluid u^. In general, also terms depending 
on the derivative of Nambu-Goldstone boson fields = d^tp, on the superfluid 
velocity and on the velocity of the director may contribute to the dissipative 
part of the stress tensor. Here the director is given by the vector pointing in the 
preferred direction. However these terms do not contribute to the off-diagonal 
components of the energy-momentum tensor for the following reasons: (1) a 
shear viscosity due to the superfluid velocity leads to a non-positive divergence 
of the entropy current [30, 31] and (2) no rank two tensor can be formed out of 
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degrees of freedom of the director if the gradients of the director vanish [15]. In 
our case the second argument is fulfilled since the condensate is homogenous and 
the fluctuations depend only on time. These degrees of freedom will generate 
additional transport coefficients, but they do not change the shear viscosities 
which we study in this letter. Thus we can write Kubo formulae which determine 
the shear viscosities in terms of the stress energy correlation functions as in (1). 
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